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ABSTRACT 



Context. Umbral Dots (UDs) are thought to be manifestations of magnetoconvection in sunspot umbrae. Recent advances in their 
theoretical description point to the need for a thorough study of their properties and evolution based on data with the highest currently 
achievable resolution. 

Aims. Our UD analysis aims to provide parameters such as lifetimes, diameters, horizontal velocities, and peak intensities, as well as 
the evolution of selected parameters. 

Methods. We present a 106-minute TiO (705.7 nm) time series of high spatial and temporal resolution that contains thousands of 
UDs in the umbra of a mature sunspot in the active region NOAA 10667 at = 0.95. The data were acquired with the 1-m Swedish 
Solar Telescope (SST) on La Palma. With the help of a multilevel tracking (MLT) algorithm the sizes, brightnesses, and trajectories 
of 12836 umbral dots were found and extensively analyzed. The MLT allows UDs with very low contrast to be reliably identified. 
Results. Inside the umbra we determine a UD filling factor of 11 %. The histogram of UD lifetimes is monotonic, i.e. a UD does 
not have a typical lifetime. Three quarters of the UDs lived for less than 150 s and showed no or little motion. The histogram 
of the UD diameters exhibits a maximum at 225 km, i.e. most of the UDs are spatially resolved. UDs display a typical horizontal 
velocity of 420 m s~^ and a typical peak intensity of 51 % of the mean intensity of the quiet photosphere, making them on average 20 % 
brighter than the local umbral background. Almost all mobile UDs (large birth-death distance) were born close to the umbra-penumbra 
boundary, move towards the umbral center, and are brighter than average. Notably bright and mobile UDs were also observed along 
a prominent UD chain, both ends of which are located at the umbra-penumbra boundary. Their motion started primarily at either of 
the ends of the chain, continued along the chain, and ended near the chain's center. We observed the splitting and merging of UDs 
and the temporal succession of both. For the first time the evolution of brightness, size, and horizontal speed of a typical UD could be 
determined in a statistically significant way. Considerable diff'erences between the evolution of central and peripheral UDs are found, 
which point to a diff'erence in origin. 

Key words. Sun: photosphere - Sun: sunspots - techniques: photometric 



1. Introduction 

The investigation of the complex fine structure of umbrae and 
] penumbrae is crucial to understanding the subsurface energy 
■ transport in sunspots. The energy transport from the solar inte- 
' rior to the solar surface outside magnetic features is mainly de- 
termined by convection, visible as granulation in images of the 
quiet photosphere. The strong and nearly vertical umbral mag- 
netic field suppresses normal overturning convection inside the 
umbra. However, it is believed that some form of residual mag- 
netoconvection is responsible for much of the remaining energy 
transport and manifests itself in the form of fine structures, such 
as light bridges (LBs) or umbral dots (UDs). In the present paper 
we consider UDs, which contribute up to 37 % of the radiativ e 
umbral flux according to lAdjabshirzadeh & Koutchmyl (Il983l) . 
DiflTerent niodels h ave been proposed to explain the umbral dots. 
IChoudhuryl (1 1 9 8 6h postulated that UDs are thin columns of field- 
free hot gas between the cluster of small magnetic flux tu bes that 
form t he subsurface structure of a sunspot according to iParkej 
According to this model, a UD is formed when an up- 
welling brings hot material in to the photo s phere. An alternative 
model has been proposed by IWeiss et aP (Il990h who consider 
UDs to be spatially modulated oscillations in a strong magnetic 
field. 

A more recen t , pro mising approach is presented by 
ISchiissler & Vogle3 (l2006l) . who used numerical simulations of 



three-dimensional radiative magnetoconvection to improve the 
physical understanding of the umbral fine structure. The simu- 
lations exhibit the emergence of small-scale upflow plumes that 
start ofif like oscillatory convection columns below the solar sur- 
face but turn into narrow overturning cells driven by the strong 
radiative cooling around optical depth unity. Most of those UDs 
show a central dark lane. The presence of dark lanes in penum- 
bral and umbr al fine structures has al r eady been observed sev - 
eral times, cf. IScharmer et al.l (l2002b : iLanghans elaP (l2007l) : 
IScharmer et al.l (l2007l). The ver ification of t he predicted dar k 
lanes in large UDs bv lBharti et a l. (2007) and lRimmeli (l2008l) . 
as well as the verification of the pre dicted photospheric strati - 
fication of bright peripheral UDs by iRiethmiiller et al.l (l2008l) . 
support the Schiissler & Vogler model of UDs. There is now a 
need to learn more about this phenomenon, with a statistically 
robust analysis of UD properties and evolution being a promis- 
ing means of achieving this aim. 

The most detailed analyses of UDs are more than 10 years 
old (Sobotkaet al., 1997a:b) and are based on data observed 
with the 50-cm SVST (Swedish Vacuum Sola r Telescope), cf . 
the recent reviews of umbral fine structures bv lSolankil(l2003h : 
iThomas & WeissI (|2004 andHo botka (2006) . The more recent 



papers of iTritschler & Schmidt (2002); .Hartkorn & Rimmel3 
(20^ andl^botka & Hanslmeier ( 2005) have concentrated on 
individual properties and lack, e.g., the determination of UD tra- 
jectories. Furthermore, the possibility of improving the spatial 
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resolution with the h elp of modern image reconst ruction algo- 
rithms is only used by lTritschler & Schmidtl (l2002h . The present 
paper aims to overcome these shortcomings, by employing data 
from the 1-m SST (Swedish Solar Telescope) equipped with an 
adaptive optics system, by restoring the data employing MFBD 
(multi-frame blind deconvolution), and determining the evolu- 
tion of UD parameters whenever possible. 

2. Observations and data reduction 

The data employed here were acquired on September 7, 2004 
with the Swedish Solar Telescope at the Observatorio del Roque 
de los Muchachos on La Palma, Spain. Technical details of the 
SST are described by Scharmer et al. (2003a). Wavefront aber- 
rations caused by the telescope and by the turbulent atmosphere 
of the Earth were partially corrected by th e adaptive optics sys- 
tem, explained in IScharmer et al.l (l2003bl) . The science camera 
was a Kodak Megaplus CCD with a pixel size of 9 yum and a 
plate scale of 0.041'' (30 km on the Sun) per pixel. The camera 
was equipped with an interference filter at the wavelength of the 
705.7 nm of the titanium oxide band head, the FWHM of this 
filter was 0.71 nm. The theoretical diff'raction limit of the tele- 
scope at the TiO wavelength is 0.18'' (130km). Due to the high 
sensitivity to umbral temperatures of TiO lines, the TiO band 
head is a good diagnostic wavele ngth range for imaging umbral 
features (iBerdyugina et al.L [20031) . A wavelength in the red was 
chosen also in order to ensure a more homogeneous time series 
due to the more benign seeing at these wavelengths. Acquisition 
lasted from 08:27 UT to 10:17 UT, i.e. a total of 110 min. The 
images were obtained in a frame selection mode that saved only 
the 8 best images of a 20- second-interval. The exposure time 
was 10 ms. The telescope pointed to the sunspot of the active 
region NOAA 10667 at cos 6 = 0.95, i.e. relatively close to the 
solar disk center (0 is the heliocentric angle). 

The data were dark current and flat field corre cted, recon- 
structed via the MFBD technique (iLofdahlL [2003b . derotated, 
destretched (iNovembeil [T9 8 8a) . and subsonic filte red with a cut- 
off' phase velocity of 5 kms"^ ([Title et al.Lll989l) . The obtained 
time series consists of 3 10 images with a spatial resolution in the 
range of ~0.18"-0.25", as we estimated from radially averaged 
power spectra. The field of view (FOV) is 37" x 59" and con- 
tains the entire considered sunspot whose umbra is divided into 
two parts by a light bridge. UDs in both parts of the umbra are 
analyzed in the next section. 

3. Data analysis 

The detailed analysis of 310 images requires an automated algo- 
rithm for the identification of the thousands of umbral dots they 
contain. A specific algorithmic challenge is the fact that UDs as 
well as the local umbral background between them cover a broad 
range of intensities. At a normal contrast (left panel of Fig. [T]) 
UDs are mainly visible near the penumbra. By displaying the 
square root of the umbral brightness instead of the brightness it- 
self numerous UDs within the dark umbral background become 
visible as well (right panel). 

Our automated UD analysis starts with isolating the umbra, 
which is done by thresholding a lowpass filtered image (averag- 
ing a squared environment of 1 1 x 1 1 pixels) at 35 % of the mean 
intensity of the quiet photosphere (Iph)- From the resulting set of 
contours we select only those longer than 3 Mm in order to avoid 
larger UDs from being connected to the umbral boundary. In the 
particular case of the studied penumbra the results are almost 
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Fig. 1. Frame taken at 09:53:20 UT: Best quality image plotted 
at normal contrast (left panel) and with increased umbral bright- 
ness (right panel). The black contour lines outline the umbra (see 
text for details). The direction to solar disk center (DC) is indi- 
cated by the arrow. Locations of umbral dots (UDs) that are dis- 
cussed further in the text are marked by the letters A-D. Minor 
tick marks are given in Mm. 



identical to identifying the two longest contours with the umbral 
boundary. The thus obtained umbral boundary is visible in the 
right panel of Fig. [T] as the black contour line. This method for 
isolating the umbral boundary automatically ensures that a local 
brightening at the end of a penumbral fibril is only considered to 
be a UD if it is isolated from the penumbra in the sense that the 
intensity between the UD and the penumbral fibril falls below 
the applied threshold of 0.35 Iph. 

In the next step the UDs in each of the 310 images (recorded 
at a cadence of 20.57 s) are detected. For this purpose several al- 
gorithms were tested, e.g. a method where, starting from the UD 
center, 8 equally distributed rays are followed until they reach 
the UD boundary which is defined as the position where the in- 
tensity drops below 50 % of the maximum intensity above the 
local umbral background. The resulting 8 boundary points lead 
to a polygon that is a good approximation of the UD bound- 
ary. The method does not work properly for UDs with a partly 
concave boundary and it cannot easily separate UDs that are 
cl ose to each other. Finall y, the multilevel tracking algorithm 
of [Bovelet & WiehJ ([200l[) . which provided the best results in 
detecting the UD boundaries, was chosen. First the MLT al- 
gorithm determines the global extrema of the umbral intensi- 
ties and subdivides this range into equidistant levels. Bovelet 
and Wiehr used MLT to distinguish between granules and in- 
tergranular lanes of the quiet Sun and found that three MLT lev- 
els are sufficient for their purpose. Since umbral dots cover a 
broad range of intensities we have to use a noticeably higher 
number of levels. We normalized our best quality image to Iph 
and found an umbral intensity range from 0.36-0.96 Iph (Note: 
This range is only valid for the best quality image, other im- 
ages may reach lower or higher umbral intensities). We found 
that 25 MLT levels is the optimal compromise between detect- 
ing as many UDs as possible and avoiding the misinterpretation 
of noise as UDs. Whereas small umbral dots, obtained with this 
choice of levels, have a typical contrast relative to the local back- 
ground of about 0.05 Iph, the noise level is about 0.005 Iph (see 
Fig. [2] for a typical intensity proffle). Starting with the highest 
intensity level all pixels are found whose intensity exceeds this 
level. This leads to several bounded two-dimensional structures. 
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Fig. 2. Illustration of multilevel tracking algorithm with 4 levels 
applied to a typical intensity profile. 



that are tagged in a unique way, which is indicated by difi'erent 
colors in the one-dimensional illustration given in Fig. [J] The 
obtained closed structures are extended pixel by pixel as long 
as the intensity is greater than the next lower level. After that 
the algorithm searches through the whole umbra again to find all 
pixels whose intensity is greater than the next lower level, which 
may lead to some newly detected closed structures. This proce- 
dure is repeated until the minimal intensity level is reached. At 
the end every umbral pixel belongs to exactly one closed struc- 
ture. The mode of operation of the MLT algorithm is illustrated 
for one dimension in Fig. [2] for the case of 4 levels. 

The minimal (iMin) and maximal (Ipeak) intensity of each 
closed structure is determined and all pixels that have an in- 
tensity lower than 50 % of this min-max range (white arrows 
in Fig. [3^) are cut. This leads to a first estimate of the UD 
boundaries (see Fig.[3]3) that are used to determine the local um- 
bral background intensities (hg), i.e. the intensities that would 
be obse rved in the absence of all UD s. We applied the method 
used by lSobotka & Hanslmeiej (l2005l) that approximates the lo- 
cal umbral background by a 2D surface fitted to the grid of 
local intensity minima , using the method of thin-plate splines 
(iBarrodale et al.L Il993h . Since the local umbral background in- 
tensities are known now (dashed line in Fig. [3]:), we determine 
the exact UD boundaries by cutting all pixels lower than 50 % 
of the maximum intensity above the local umbral background 
(see white arrows in Fig. [3]:). Fig. [3]1 illustrates that the result- 
ing UD boundaries are similar to our first estimate in Fig. [3t 
so that we don't need further iterations (this was the case with 
most identified UDs). Employing this procedure we found, on 
average, 323 UDs per image. Sometimes, our algorithm recog- 
nizes strong elongated bright structures as UDs that we would 
not consider as a UD by visual inspection of the images. Via a 
spot-check on selected images we estimate that the number of 
misidentifications is lower than 1 %. Due to the large number of 
detected UDs we are not able to remove the misidentifications 
by hand and accepted them as noise. 

Since a UD is an extended structure we determine the coor- 
dinates of the brightest pixel (peak intensity) and save them as 
the UD's position. This method is applicable because the noise 
was sufl&ciently reduced by our subsonic filter, as demonstrated 
by a typical intensity profile in Fig. [2l We also determine the 
UD's diameter, defined as the diameter of a circle of area equal 
to that within the boundary of the UD. After the positions of all 
UDs of every image are known, the motions and trajectories of 
UDs are determined. The continuation of a trajectory in the im- 
age at the next (previous) time step is determined by finding the 
UD that is closest to the UD's current position. If no UD can 
be found within a 5 pixel neighborhood (theoretical diff'raction 
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Fig. 3. Determination of the boundary of umbral dots. 

limit) of the current UD then the tracking stops. In a loop over all 
images every UD is tracked backward in time until its birth and 
forward in time until its death. The tracking must be tolerant to 
the occasional image with lower image quality in which the UD 
may not be correctly identified (specially the smaller and fainter 
ones). In practice we allow for a gap of up to two images. If a 
UD is present at nearly the same location on both sides of the 
gap, the tracking is continued. In this manner we found 12836 
UD trajectories that are weakly lowpass filtered in space (aver- 
aging the positions of the same UD in 15 consecutive images) 
in order to further reduce seeing-induced noise. From now, we 
call such a smoothed trajectory simply tra jectory. We note that 
5949 of the 12836 UDs are only identified in a single image. We 
decided not to ignore them, because these bright dots in the um- 
bra are detected rather well, even if only for a very short time. 
We assigned a zero trajectory length and a lifetime of 20.57 s to 
these 5949 UDs. 

4. Results 

4.1. Qualitative results 

A first impression of the temporal evolution of the smallest um- 
bral structures is reached by making a movie of the reconstructed 
time series of images. Some interesting phenomena are found by 
the visual inspection of this movie and are explained briefly be- 
low. 

The sunspot has two umbrae, a smaller and a roughly twice 
larger one separated by a light bridge (LB). The LB contains 
a clearl y visible dark lane in agreem ent with the observations 
of, e.g., iBerger & Berdyuginal (l2003l) . The lane is closer to the 
limb ward edge of the LB. Possibly this is because the sunspot 
was observed at a heliocentric angle of ^ = 18°, so that projec- 
tion e fl^ects may cause the observed asymmetry fe.g. lLites et al.L 
|2004|) . However, the LB also displays another major asymme- 
try: the movie exhibits many UDs that are born within the LB 
and move into the larger umbra, i.e. towards the solar disk cen- 
ter (arrow in Fig. [T]), while almost no UD leaves the LB into 
the smaller umbra, i.e. towards the solar limb (antiparallel to the 
arrow). 

The data clearly show that bright UDs often form chains. The 
most prominent chains often start from a penumbral filament and 
UDs are found to move along the chain until they dissolve. The 
horizontal motion of the UDs is preferentially along the chain 
and is directed from the ends of the chain to its center, where 
the UDs disappear. Fig. |4] displays the most prominent chain of 
UDs in the observed umbra. The length of this chain is about 
3200 km and it is about 350 km wide. In the left panel, from 
08:58:56 UT, the chain appears as a simple succession of UDs. 
One can see the same region in the right panel 28 minutes later. 
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Fig. 4. Most prominent chain of umbral dots in the studied um- 
bra, located near label A in Fig. [T] (FOV is 5.4 x 7.4 Mm). 
The left panel shows the chain composed of several typical UDs 
along a curved line at 08:58:56 UT. The right panel shows the 
same chain 28 min later. At this later time the lower half of the 
chain looks similar to a narrow light bridge. The typical dark 
lane is clearly recognizable. 



08:53:22 


08:56:44 


09:00:09 


09:03:18 












% 







Fig. 5. Splitting of an umbral dot near the position B in Fig. [T] 
(FOV is 1.2 X 1.2 Mm). The red and yellow lines are the UD 
boundaries as detected by the method explained in the main text. 



The appearance of the chain changed and now the lower part 
of the chain looks similar to a narrow light bridge. The typical 
dark lane of a LB can be seen clearly, even if we degrade the 
image quality to the lower level of the left panel by convoluting 
the image with a point spread function of a circular pupil (not 
shown). This degradation was carried out to compensate for the 
higher spatial resolution of the later image. High spatial resolu- 
tion is demon strated by the presence of dark-cored penumbral 
filaments (see iLanghans et al.l l2007l) . which are generally ob- 
served only at the highest resolution in the blue. The presence of 
the dark lane is an indication that the chain eventually evolved 
into a fu lly developed LB so me hours after the end of the record- 
ing dKatsu kaw a et a l. j2007]) . 

In some exceptional cases we observe the splitting of a single 
UD into two parts that continue their life as independent UDs. 
Such a splitting is displayed in Fig. [51 It concerns a UD located 
close to the penumbra, near point B in Fig.[T] The opposite case, 
the merging of two UDs into a single one, can also be observed 
in a few rare cases: See Fig. [6] for an example. Nearly identical 
phenomena in the temporal evolution of penumbral grains were 
observed by Hirzberger et al. ( 2002). 

An interesting sequence of events is illustrated in Fig.[7l First 
the merging of two UDs can be observed followed immediately 
afterwards by the splitting of the resulting, unified UD into two 
new UDs. The sequence looks similar to an elastic impact in 
a Newton pendulum consisting of two balls: The nearly motion- 
less UD 1 is hit by the moving UD 2. The impact of the two UDs 
brings the second UD to a standstill while the first UD starts to 




Fig. 6. Merging of two umbral dots near position C in Fig. [T] 
(FOV is 1.2 X 1.2 Mm). 




Fig. 7. Elastic collision of two umbral dots. The sequence of 
events is seen near the position marked D in Fig. [T] (FOV is 
0.9 X 0.9 Mm). UD 2 moves towards UD 1 and merges with 
it, before separating from it again. 



move roughly in the direction of the first UD. Note that this is 
simply an empirical description and we do not propose that this 
is what physically happens (or that the unified UD breaks up into 
the same parcels of gas which united to form it). 

If two UDs come close to each other then the visual impres- 
sion of a single UD exhibiting a dark lane can occur, see last 
panel of Fig. [7] for an example. None of the UDs in our data 
set seems to stay in such a state for a significant fraction of the 
UD lifetime. Since t he dar k lanes as seen in the simulations of 
ISchiissler & Vogled (l2006h are visible for most of the UD life- 
time, we conclude that we find no clea r evidence for such dark 
lanes. Th i s diff'e rence to the results of iBharti et al.l (l2007l) and 
iRimmelj (l2008l) may stem from the diff'erent wavelengths of the 
analyzed data. The wavelengt h can have a remarkabl e eff'ect on 
the detected fine structure (e.g. lZakharov et al lHooi). 

The properties of a UD change along its trajectory. In order 
to assign a property to the entire trajectory we either average over 
all points along the trajectory (which is expressed by introduc- 
ing an upper index "Mean") or we determine the maximum value 
reached by that parameter over all trajectory points (which is ex- 
pressed by the upper index "Max"), e.g., D^^^^ means the UD 
diameter averaged along the trajectory, while /^^^ is the maxi- 
mum value of all peak intensities along a UD trajectory (where 
the peak intensity is the largest intensity within the UD boundary 
at a given point in time). 

Fig. [8] shows the best image with all identified UDs marked 
by circles. The circles are centered on the positions of i.e. 
the position of maximum intensity. In the left panel the circles 
have a constant radius, while their radii are proportional to /^^^ 
in the right panel. As a result we get an impression of the spatial 
distribution of UD occurrence as well as of the spatial distri- 
bution of UD brightness. Obviously, there is hardly any part of 
the umbra which does not support umbral dots. Only very local- 
ized small voids are visible. The brightness distribution of UDs, 
however, is rather inhomogeneous, with clear concentrations of 
bright UDs and regions harboring mainly dark UDs (mainly in 
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Fig. 8. Spatial distribution of UD occurrence (left panel) and UD 
maximum brightness (right panel). Each UD is plotted only once 
along its entire trajectory, at the position of its maximum bright- 
ness. 
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Fig. 9. Dependence of UD filling factor (dotted blue line) and 
granulation contrast (solid black line) on the interval of time, t, 
since the first image of the series was recorded. 

the upper right part of the upper umbra and in the lower part of 
the lower one). Since the UDs of a chain (like the chain close to 
label A in Fig. [T] and shown in more detail in Fig. [4]) appear to 
move along the chain and hardly in the direction perpendicular 
to it, there is often a narrow void directly beside the chain. 

4.2. Quantitative properties 

The filling factor, i.e. the sum of all UD areas relative to the total 
umbral area, is correlated to the image quality. Essentially, the 
filling factor is constant over the entire period of observations 
(see Fig. [9]), which is important for the later determination of 
the time dependence of their properties. On average we deter- 
mine a value of about 11%. Fig. [9] shows also the image contrast 
from an undisturbed granulation area outside the sunspot which 
demonstrates the high homogeneity of the image quality in our 
time series. 

The histogram of the UD lifetimes is displayed in Fig. [TOl 
with a logarithmic y axis. One can see that most of the UDs 
live for a short time. These short-lived UDs move over short dis- 
tances which leads to physically nonsensical velocities due to 
the discretization of lifetime and distance. Consequently, when 
discussing trajectories and velocities of UDs we only consider 
the 2899 trajectories of UDs with lifetimes greater than 150 s. 
The histogram is nearly linear for lifetimes between 5 and 60 
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Fig. 10. Histogram of lifetimes T of all 12836 UD trajectories 
(solid line) and their mean and median value (dashed lines and 
text labels). The bin size is 300 s. 



minutes, which, due to the logarithmic vertical scale, suggests 
an exponential distribution of lifetimes. The excess of UDs with 
short lifetimes may partly be due to seeing. Also given in Fig.fTOl 
is the mean lifetime (180 s) and the median lifetime (41 s). If we 
consider only the 2899 trajectories of UDs with lifetimes greater 
than 150 s then we find a mean lifetime of 630 s and a median 
of 390 s. Note that 281 UDs are already present in the first im- 
age and 344 UDs are still present in the last image, whereas only 
one UD survives the whole sequence. Ignoring those 625 UDs 
reduces the mean lifetime from 180 s to 152 s. The median life- 
time as well as the shape of the histogram do not change, because 
the number of incomplete UD trajectories is small compared to 
the total number of UDs. Thus we decided to neglect this eff'ect 
and consider all 12836 UDs in the following text. 

The histogram of the mean UD diameters is plotted in 
Fig.[TT]a. The mean UD diameters vary between 50 and 750 km. 
The UDs have a mean diameter of around 229 km and 95 % of 
the UDs are spatially resolved at our diff'raction limit of 130 km. 
ISobotka & Hanslmeied (l2005h found 175 km for the mean UD 
diameter. They determined the UD boundaries by finding all 
pixels with downward concavity, whereas we used all pixels 
whose intensity is greater than 50 % of the Ipeak - hg range, 
^hhresh = (Ithresh " hg)l{heak " hg) = 0.5. As onc Can see in 
Fig. [12] we would also find a mean UD diameter of 175 km if a 
threshold of about Althresh = 0.67 would be used. Nevertheless, 
we prefer to continue our analysis with a value of 0.5 in anal- 
ogy to the FWHM definition. Clearly, the employed thresh- 
old influences the filling factor as well, roughly quadratically. 
Remarkably, the histogram is nearly symmetric, which supports 
the conclusion that most UDs have been resolved. 

The histogram of the mean horizontal velocities, i.e. the quo- 
tient of trajectory length and lifetime, is plotted in Fig. [TT] b 
and exhibits a broad distribution from to more than 1 km s"^ 
with a significant maximum at 350 m s"^ . The velocity distribu- 
tion is slightly asymmetric with a small tail to higher velocities. 
Fig.[TT]c shows the histogram of the mean peak intensities, i.e. 
the mean of all peak intensities of the points along the trajectory. 
All intensities are normalized to the mean intensity of the quiet 
photosphere (Iph)- Just 3 of the 12836 UDs reach a brightness 
greater than that of the quiet Sun, whereas most of the UDs are 
about half as bright as the quiet photosphere. The distribution 
is asymmetric, with a tail to higher intensities. Note that these 
brightnesses are strongly wavelength dependent and can not be 
easily compared with values published in the literature ( Sol ankiL 
l2003l) . Lastly, the histogram of the distances between the UD's 
birth and death position (Lbd) is given in Fig.[TT]d. With increas- 
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Fig. 12. Mean UD diameter averaged over all trajectories as a 
function of the intensity threshold Alfhresh (defined in main text) 
that is used to determine the UD boundaries (see also Fig. [3] and 
its explanation). 
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Fig. 13. Scatterplot of mean UD diameter versus mean peak in- 
tensity. The solid green line connects binned values. 



ing distance the number of UDs decreases exponentially, so that 
only a few UDs travel over long distances in their life, most UDs 
do not leave the vicinity of their birth position. The maximum 
observed birth-death distance is 2 Mm, about 20 % of the upper 
umbra's diameter of roughly 10 Mm. 

A UD size versus UD peak intensity scatterplot (Fig. [T3]) 
reveals that there is a weak correlation between UD size and 
brightness. The brightest UDs are not the biggest ones and 
large UDs are not the brightest ones. The solid green line con- 
nects binned values (obtained by averaging 100 data points 
with similar intensities) and shows that bright UDs are on 
average a bit lar ger than small ones. The relation found by 
iTritschler & Schm idt ( 2002) is qualitatively confirmed, although 
they only considered UD intensities in individual snapshots 
while we tracked the temporal development of the UDs over 
their lifetimes. 

The relation between a UD's mean size (i.e. the size averaged 
over the lifetime) and its lifetime is plotted in Fig. [14] (the green 
curve is obtained after binning over 100 data points). The binned 
values show an increase in D^^^^ with increasing lifetimes for 
short-lived UDs. For the longer lived ones size and lifetime do 
not correlate. The UD sizes scatter more for short lifetimes. All 
long-lived UDs are of intermediate size of around 290 km. The 
large, short-lived UDs are all present in the first image of the 
time series, so that their lifetime would actually be larger if we 
had started our observation earlier. 

In the li terature we often find a separa tion into two UD 
classes, e.g. iGrossmann-Doerth et al.l (Il986h find a diff'erence 



Fig. 14. Scatterplot of mean UD size versus lifetime. The solid 
green line connects binned values. 



between peripheral and central UDs, i.e. between UDs that are 
born close to the umbra-penumbr a boundary and UDs that are 
bor n deep in the umbra, whereas iHartkorn & Rimmel3 (|2003|) 
an d ISobotka et all ([T997bh distinguish bright and dark UDs. We 
use diff'erent properties to find reasonable distinctions between 
types of UDs. E.g. from now on a given UD is called a periph- 
eral UD (PUD) if the UD's birth position is closer than 400 km to 
the umbral boundary, otherwise it is termed a central UD (CUD). 
(The selected threshold comes from the histogram of the dis- 
tances between the UD's birth position and the umbral bound- 
ary (not shown) which shows a maximum at around 400 km.) 
Alternatively, if the distance traveled between birth and death 
position is larger than 750 km then we call it a mobile UD, oth- 
erwise a stationary UD. (We plotted the trajectories of all UDs 
whose Lbd was greater than a threshold, which was determined 
by starting at a small value and increasing it step by step. We 
stopped at 750 km which is the smallest Lbd at which no trajec- 
tories occurred anymore in the central part of the umbra.) The 
aim here is not to separate UDs into distinct classes by a single 
property, e.g. a histogram of Lbd (Fig. Eld) does not show two 
peaks, even if restricted to long-lived UDs. However, as we shall 
see below the most mobile UDs are formed near the penumbra, 
while the least mobile ones are mainly formed deep in the um- 
bra. Such a distinction may help to guide theory towards a better 
understanding of the origin and evolution of UDs with diff'erent 
properties and at diff'erent locations. 

The strong concentration of bright UDs near the umbral 
boundary (or along proto-light bridges) is already clear from 
Fig. [8] (right panel). Fig. [15] shows a separation by birth-death 
distance. The left panel displays the longest UD trajectories, i.e. 
only the mobile UDs are shown. These UDs are also relatively 
long-lived. (Smallest lifetime of this UD class, that contains 85 
UDs, is 15 min.) The right panel shows trajectories of long-lived 
UDs with a small birth-death distance. A clear separation by the 
birth position is readily identifiable. Almost all trajectories with 
a large birth-death distance start close to the umbra-penumbra 
boundary and these UDs move nearly radially into the umbra. A 
visual inspection of the movie of the reconstructed time series 
of images shows that many of these UDs are former penumbral 
grains that broke away from the penumbra. In contrast, many of 
the UDs with a small birth-death distance are born in the um- 
bral interior and move along a closed loop or jitter around their 
birth position. We cannot say if this jitter has a physical back- 
ground or if it is caused by residual seeing-induced noise. As 
mentioned in section 3 we determined the umbral boundary in- 
dividually for each image but we show only that corresponding 
to the best quality image as black contour line in Fig.[T5l[T71 and 
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Fig. 11. Histogram of mean diameters (a), mean horizontal velocities (b), mean peak intensities (c), and distances between birth and 
death position (d). (a), (c), and (d) are plotted for all 12836 UD trajectories and (b) for the 2899 trajectories of UDs that lived longer 
than 150 s. The location of the maximum, the mean, and the median of the distribution is indicated in each frame. The bin sizes are 
30 km for (a), 100 m s"^ for (b), 0.03 Iph for (c), and 90 km for (d). 



[T8l Consequently some trajectories (or parts of them) are out- 
side the black contour line, although they are always inside the 
umbra at the time of their occurrence (see, e.g., the bottom-left 
corner of the left plot of Fig.fTSl). 

Let us now consider more quantitatively the fact that the 
mobile UDs prefer to move radially towards the umbral cen- 
ter. Fig. [16] displays a histogram of the UD's deflection angles 
adefl that is defined as angle between the line connecting the um- 
bral center and the UD's birth position and the line connecting 
the UD's birth and death position. Radially directed inward flow 
will lead to adefi = 0° and an outward flow to Ofdefl = 180°. 
Obviously, this definition makes sense only for the 5689 UD tra- 
jectories that have diff'erent birth and death positions. The solid 
black line in Fig. [161 shows the histogram for all these UDs and 
exhibits a clear tendency for a radially directed inward motion 
(42 % of the UDs are found to have a deflection angle lower than 
45°). This tendency is much more significant if we only consider 
mobile UDs, see the dotted blue line (88 % of the mobile UDs 
are found to have a deflection angle lower than 45°). However, 
histograms calculated for central UDs and for peripheral UDs 
(not shown) lead, in principle, to the same shape as for all UDs, 
i.e. UDs born close to the penumbra do not show a significantly 
higher tendency of radially inward directed motion than the cen- 
tral UDs. 

The trajectories of large UDs are drawn in Fig. [T7] (left 
panel). They are born throughout the umbra, with a tendency 
to cluster in the darker part of the umbra. The vicinity of the 
light bridge is avoided. This UD class contains long trajectories 
as well as short ones. The trajectories of the brightest UDs can 
be seen in the right panel of Fig.[T71 These UDs all emerge close 




ttdefi [deg] 

Fig. 16. Histogram of radial deflection angles Ofdefi of all UD 
trajectories (solid black line) and of the mobile UD trajectories 
(dotted blue line). 



to the penumbra, the light bridge, or the prominent UD chain 
(label A in Fig.[B. 

Many UDs with a preferred direction of motion arise near 
the light bridge. Most of the UDs that emerge on the disk cen- 
ter side of the light bridge (i.e. into the large umbra) move away 
from the light bridge while many of the limb side UDs (i.e. those 
formed in the small umbra) move towards the light bridge (see 
Fig. [18]). A high density of UDs is formed by splitting ofl' the 
light bridge, but all in one direction, in which the LB is corru- 
gated and unsharp (disk center, large umbra side), while on its 
other straight and sharp side nearly no UDs leave the LB. 

The mean diameter D^^^^, mean horizontal velocity 
^Mean^ mean peak intensity mean intensity contrast 
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Fig. 15. UD trajectories in the upper umbra. The bright ends of the trajectories denote the positions of the UD's birth and the dark 
ends show the position at death. The black contour line corresponds to the umbral boundary as detected in the best quality image. 
The left panel shows all UDs whose distance between the birth and the death position was greater than 750 km, the right panel 
shows all UDs with a birth-death distance smaller than 150 km whose lifetime was greater than 800 s. 




Fig. 17. UD trajectories in the upper umbra. The left panel shows all UDs with mean diameter greater than 350 km and the right 
panel shows all UDs with a mean peak intensity greater than 0.65 I ph. 



{Ipeaklhg)^^^^, lifetime T, birth-death distance Lbd, trajectory 
length Ljraj (see definitions of these quantities in the previous 
text), and the number of UDs of the types or classes men- 
tioned earlier in this section are summarized in Table [T] We av- 
eraged over all trajectories of a UD class. The standard devia- 
tion cr is given after each average value. As mentioned earlier, 
we consider only UDs that lived longer than 150 s in all cases 
in which we calculate the mean velocity; only the first line in- 



cludes all UDs. The last two rows consider only UDs that are 
close to chain A or the light bridge, respectively. According to 
Table [T] the diff'erence between peripheral UDs and central UDs 
is not so large (using the simple categorization described above). 
The largest relative diff'erence is in the brightness. For all other 
parameters the diff'erence is less than Icr. In contrast to that, the 
diff'erence between mobile UDs and stationary UDs is more sig- 
nificant. A relatively large diff'erence of more than Icr is found 
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Table 1. Characteristic values of UD parameters for different UD classes defined in the main text. 
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Fig. 18. UD trajectories near the light bridge with birth-death 
distance greater than 300 km. 



for the mean horizontal velocity, for the mean intensity contrast, 
for the lifetime, and, again for the brightness. On average mo- 
bile UDs are brighter, they move faster, and they live longer than 
stationary UDs but they have similar sizes. 

Sometimes one can observe complete chains of successive 
UDs that are very close to each other (see label A in Fig. [T] 
and the left panel of Fig. 14]). According to Table [T] these UDs 
are relatively bright and they move along the chain from the 
endpoints of the chain towards its center, where they disap- 
pear. Their mean peak intensity is 0.68 /^/^ which is significantly 
higher than 0.51 Iph, the average of all UDs, but is comparable 
to that of the peripheral and mobile UDs. The mean diameter as 
well as the mean velocity of the UDs within the chain is slightly 
above average. Table [T] also reveals that UDs that are born close 
to the light bridge show on average a significantly higher bright- 
ness and contrast than the mean UD but all other properties do 
not show remarkable diff'erences. 

The umbral background is brightest near the penumbra and 
gets darker towards the center of the umbra. Fig. [19] shows that 
the UD peak intensity correlates with the umbral background in- 
tensit y, with the mean ratio Ipeaklhg = 1.2 + 0.1 (intensity con- 
trast). [Sobo^a^Hannmeie^ (llOQ^ found an intensity contrast 
of 1.8 for a wavelength of 451 nm and 1.6 for the wavelength 
602 nm. Obviously, the intensity contrast not only depends on 
the wavelengths but also on the spatial resolution that can be 
slightly diff'erent even if we compare data from the same tele- 
scope. Additionally, the intensity contrast can also be aff'ected by 
the subsonic filter and the MFBD image restoration we applied. 
Importantly, however, both binning the points in Fig. [19] (solid 
green curve) and a linear regression (dotted red curve) indicate 
that the contrast increases nearly linearly with 4^. A scatter plot 
of the UD intensity contrast versus the shortest distance between 
the UD birth position and the penumbra (not shown) reveals that 
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hg [Iph] 

Fig. 19. UD peak intensity versus umbral background intensity. 
The solid green line connects binned values. The dashed cyan 
line displays the theoretical, lower limit of the peak intensities. 
A linear fit to the data results in the dotted red line. 

this UD contrast is not constant over the umbra. The closer the 
UD is born to the penumbra the higher its intensity contrast, al- 
though the contrast does not drop as rapidly from the umbral 
boundary as the UD brightness does, so that partly, the depen- 
dence on distance is due to the dependence on I^g (Fig. [19]). We 
also found that on average the long-lived UDs have a higher con- 
trast than the short-lived ones. 

In a further step we subdivided the umbra into several boxes 
and determined the probability that a UD is born in such a box. 
The map obtained in this manner shows a uniform distribution of 
the UD birth probability (not plotted). Only for very small box 
sizes do the dark umbral nuclei (see left panel of Fig. [8]) become 
visible as locations of reduced UD production. 

Finally, we are interested in the temporal evolution of the UD 
properties over their lifetimes. To this end we normalize all UD 
lifetimes to unity and average the temporal evolution of the di- 
ameters, peak intensities, intensity contrasts (Ipeak/hg), and hor- 
izontal velocities of the 2899 trajectories with T > 150 5". In or- 
der to weight all trajectories equally, i.e. independently of their 
lifetime, we up-sample all trajectories to 310 points of time via 
interpolation. (Our time series contains 310 images, so that no 
trajectory can consist of more than 310 points.) Then we aver- 
aged the UD parameters with the help of our binning method, 
which is applied separately for the 621 PUDs and the 2278 
CUDs. Each bin contains 15000 points for PUDs and 50000 
points for CUDs. The results are plotted in Fig. [20] and show 
that the mean PUD is smaller, brighter and moves faster than the 
mean CUD (as could already be deduced from Table [T]). More 
importantly, there are distinct diff'erences in their mean evolu- 
tion. Whereas both types of UDs share the property that their 
diameters evolve rather moderately over time (the increase af- 
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ter birth and the decrease before death are less than 10 % of the 
maximum diameter (see panel a), the evolution of their bright- 
ness (panel b) and in particular of their contrast (panel c) dif- 
fer considerably. While the mean CUD displays an initial gentle 
brightening followed by an equally gentle darkening, the mean 
PUD darkens continuously. The small magnitude of the change 
in brightness may be due to the fact that we have averaged over 
UDs with very different absolute intensity. More information 
may be gleaned from the contrast, i.e. the peak intensity divided 
by the local umbral background intensity, plotted in panel (c). 
The mean PUD initially remains almost constant, exhibiting a 
slight maximum at around 1/3 of the mean lifetime before drop- 
ping rapidly over the remaining portion of its life. The contrast of 
the CUDs displays a much more symmetric evolution. The birth 
velocity of the mean PUD is nearly 50 m s"^ higher than for the 
mean CUD, while the velocity at death of the two UD types is 
similar, see panel (d). Both velocity curves show an initial in- 
crease, followed by a decrease. As in the case of the contrast the 
velocity profile is much more symmetric for the CUDs. 

5. Discussion and conclusions 

We have analyzed a time series of images of a mature sunspot 
close to solar disk center. Due to the excellent image qual- 
ity we were able to resolve thousands of U Ds. Exhaustive 
UD an al yses can be found in earlie r papers (ISobotka et aL . 
1997a[lbl: iHartkorn & RimmeleL l2003l: ISobotka & Hanslmeiei . 
20051: ISobotkaL 120061) . but the present article is the first detailed 
UD study of a long time series of reconstructed images with the 
consistent high resolution of a 1 -meter telescope. 

Trajectories, lifetimes, diameters, horizontal velocities, peak 
intensities, and distances between birth and death locations were 
determined by tracking single UDs over the time series. These 
characteristic values were used to look for reasonable separa- 
tions into UD classes. In the following we summarize the ob- 
tained results and compare them with other investigations in the 
literature: 

1 . There is hardly any part of the umbra which does not sup- 
port UDs, but the UD brightnesses depend strongly on the 
location wit hin the umbra , which confirms the previous ob- 
servation of ISobotka eTaP (Il997bl) . 

2. The histogram of lifetimes shows an exponential distribu- 
tion, i.e. a UD does not have a typical lifetime. More than 
3/4 of all studied UDs lived less than 150 s and their motion 
was negligible. The exponential distri bution is in qualitative 
agree ment with the res ults obtained by Sobotka et al.l (Il997al 
Il999h . Quantitatively, ISobotka et aP (Il997al) obtains a me- 
dian lifetime of 6 min for an umbra of about 6 Mm diameter 
and a median of 12 min for a 4 Mm pore (ISobotka et al.L 
119991) . whereas we find a median value of 0.7 min for a 
roughly 10 Mm umbra. Note that these median values de- 
pend strongly on algorithmic constraints as well as on the 
cadence of the time se ries. For example, the method used by 
ISobotka et al.l (Il997al) cannot lead to lifetimes shorter than 
1.5 min. If we only consider UDs with lifetimes greater than 
1.5 min our median increases to 4.1 min. Irrespectively of 
which of these two values we use, our results a re consis- 
tent with the conclusion of ISobotka et al. I (Il999h that UDs 
are more stable in a weak magnetic field if we assume a 
direct correlation between umbral diameter and magnetic 
field strength (cf . Kopp & Rabin, 1992). Alternatively, due to 
the strong dependenc e of umbral brightness on umbral size 
(iMathew et al.L l2007l) the lifetime may be influenced mainly 



by the radiative flux or umbral temperature. The scatterplot 
of the mean umbral background intensities versus the UD 
lifetimes (not shown) is consistent with both possible expla- 
nations mentioned above: UDs live longer in brighter parts of 
the umbra. Due to the non-linear, monotonically decreasing 
relation between mag netic field strength and background in- 
ten sity as observed by Kodd & Rabin (1992') an d confi rmed 
by iMartmez & Vazq uez (1993); Solanki et alj (Il993h . this 
implies that UDs live longer in regions of weak field. 

3. The histogram of mean diameters exhibits a maximum at 
225 km (0.3 r') and descends from there towards the diffrac- 
tion limit, so that we expect the majority of UDs to have been 
spatially resolved. This seems not to be the case in many 
of the earlier papers because there a mono tonic decrease 
was obtained towards higher diamete r s (see ISobotka et al.L 
Il997ai Il999h . ISobotka & Hanslmeied (l2005l) also analyzed 
data obtained with the 1 -meter SST. These data lead to a his- 
togram that is qualitatively similar to ours. The mean diame- 
ter of 175 km (0.24''), as well as the average filling factor of 
9 % is, however, noticeable smaller. This diff'erence can be 
explained by the use of a diff'erent method to determine the 
UD boundary. An increase of our brightness threshold to de- 
termine the UD bou ndary leads to small er diameters and to 
lower filling factors. iHamedivafal (l2008l) used an improved 
method of image segmentation and also found a mean diam- 
eter of 230 km and a similar shaped histogram. 

4. The mean horizontal velocity of those of our UDs that live 
longer than 150 s is 420 m s"^ which is significantly higher 
than the 210 ms"^ reported by iMolownv-HorasI ([1994) and 
higher than the 320 m s"^ found by ISobotka et al.l (1 19991) . In 
both studies the mean horizontal velocity was calculated by 
means of least- squares linear fits of the x and y coordinates 
of all trajectory points, which leads to an underestimation of 
velocity in case of curved trajectories. Our histogram of hor- 
izontal velocities shows a maximum at 350 ms"^ whereas 
some UDs can reach velocities above 1 kms ~^. This is i n 
qualitative agree n ient w ith th e histograms of iKitail (Il986l) . 
iMolownv-HorasI (Il994l) . and 'Hamedivafa" ('2008') but dis- 
agrees with the results o fjobotka et al. (1997b, 1999) whose 
histograms do not exhibit a maximum; they peak at zero ve- 
locity and show a monotonic decrease toward higher veloci- 
ties of up to 1 km s"^ . The majority of our UDs moves irregu- 
larly around the birth position. However, there are some mo- 
bile UDs that travel over long distances within their lifetime. 
Almost all mobile UDs emerge close to the umbral border, 
i.e. near the penumbra, they are brighter than the average and 
their horizontal motion is preferentially directed towards the 
center of the umbra. The mean velocity of our mobile UDs is 
680 m s~ | , which is in good agree ment with the recent obser- 
vation of lKatsukawa et al.l (l2007l) who found a mean velocity 
of 700 m s-^ 

5. The relation between mean UD size and lifetime is non- 
linear. On average, the si ze of UDs inc reases with lifetime, 
which was also found by ISobotka et al.l ([T997a). but in con- 
trast to their work we find a narrow size distribution of 
around 290 km for long-lived UDs. 

6. UDs that were born close to the penumbra show a 
significantly higher contrast than the UDs of the um- 
bral interior. The mean UD intensity contrast Ipeak/hg is 
1.2 which is smaller than t he value of 1.6 reported by 
ISobotka & Hanslmeiei (l2005h . This may partly be due to the 
longer wavelength of our observation. Additionally, our sta- 
tistical ensemble contains many more UDs. In particular we 
took many UDs with low contrast into account, made pos- 
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Fig. 20. Temporal evolution of the UD diameter (a), UD peak intensity (b), intensity contrast (c), and velocity (d) separated for 
peripheral UDs (solid red line) and central UDs (dashed blue line). The UD lifetime is normalized to unity. 



sible by the multilevel tracking technique (iBovelet & Wiehii 
I2OOII) we employed to identify UDs. Consequently, we be- 
lieve that the lower UD contrast we find is not due to a 
lower resolution or higher stray light, but rather to diff'er- 
ences in identification of UDs and in particular the diff'er- 
ence in wavelength. We stress that the UD contrast itself 
depends on the background intensity; the higher the inten- 
sity the stronger the contrast. Also it cannot be ruled out that 
there could be systematically diff'erent contrasts between dif- 
ferent sunspots due to intrinsicly diff'erent physical proper- 
ties of the spots. 
7. Whereas the temporal variation of the UD diameter is 
qualitatively similar for UDs formed close to the penum- 
bra (PUDs) and those formed in the body of the umbra 
(CUDs), their intensity contrast and horizontal velocity dis- 
play contrasting evolutions. The mean PUD shows a con- 
tinuous darkening which is in agreement with the results of 
lKitaietal l(l2007^for a single typical PUD. The typical CUD 
of lKitaietal.1 (l2007h is found to increase in brightness lin- 
early and then to darken linearly with time, whereas our re- 
sults for the mean CUD show a non-linear increase in bright- 
ness until nearby half of the lifetime followed by an again 
non-linear decrease. The clear diff'erence between PUDs and 
CUDs in the behavior of their contrast and mean velocity 
may b e a result of the di ff'erent origin of the two types of fea- 
tures (iKitai et al.l l2007l) . We confirm from visual inspection 
of a subset of PUDs that PUDs are formed when penumbral 
grains cross the umbral boundary. 

A comparison of the results with the simulations of 
ISchtissler & Vogle3 (l2006h shows a better agreement with 
CUDs, than PUDs. For example, the simulated UDs display a 
gradual increase in contrast followed by a gradual decrease, just 
as CUDs. They also display little proper motion. This quali- 



tative agreement further strengthens the interpretation of UDs 
as localized columns of o verturning convection proposed by 
ISchiissler & Vogled (l2006h . Hinode dat a had earlier sugge sted 
the presence of dark lanes in large UDs (iBharti et al.L l2007l) and 
revealed a decrease in the magnetic field strength with depth, as 
well as a n upflow associated with a temperature enhancement 
(Riethm iiller et al.l l2008l) . in good qualitative agreement with 
the simulations. A detailed analysis of the simulations similar 
to the one carried out here would allow a more quantitative com- 
parison. 

The PUDs have significantly diff'erent evolution histories 
than the simulated features. They start at a higher speed (iKitaiL 
and in particular display their maximum br i ghtne ss right 
after the beginning of their fife (cf. iKitai et al.L l2007h . This, 
combined with the fact that they are born very close to the 
penumbra, or actually b y breaking away from the penumbra 

omas &Weissl.l2054h . and move radially towards the umbral 
center (IKitaiL I 119861) supports that these are two distinct types of 
UDs based on their origin and evolution, although their physical 
structure is relatively similar (iRiethmliller et al...,2QQ8.) . 
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